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We consider multi-Galileon theory, the most general Galilean invariant theory with N scalar fields 
linearly coupled to the trace of the stress-energy tensor. We study the behavior of perturbations on 
' a static spherically symmetric background with a massive point source, and show that, under the 

assumptions of stability and successful Vainshtein screening, solutions cannot be found that are free 
of both superluminal propagation and slowly moving, strongly coupled fluctuations. The latter imply 
that the theoretical and phenomenological issues related to a very low strong-interaction scale cannot 
be avoided in this model. 



I. INTRODUCTION 

Galileon theory was introduced as a local infrared modification of general relativity (GR) that could 
potentially explain the observed cosmic acceleration in a natural way, thanks to the existence of stable 

m : 

, self-accelerating background solutions QJJ] . Its origin can be traced back to the DGP model 120, a higher- 

dimensional brane-world model in which the graviton appearing in the 4-dimensional effective theory prop- 
C\l . agates additional degrees of freedom. It is possible to thoroughly study most of the interesting properties 

^ ' of this setup in the so-called decoupling limit |3|], in which the modifications to GR are encoded in a scalar 

field 7T whose Lagrangian is invariant under the Galilean symmetry -k — > -k + b^x^ + c. This motivated 
the study in Ref. QJJ] of the most general theory exhibiting this Galilean invariance I40 (with second-order 
equations of motion), providing a framework in which well-behaved self-accelerating solutions could be 
found U]. 

The phenomenology of Galileon theories has proved to be extremely rich, and in the lapse of only a 
few years this class of theories has come to be recognized as an important alternative to other more usual 
modified gravity theories. Galilean invariance has found applications in a variety of models relevant to 
cosmology [(]] and astrophysics [7]. In particular, the fact that Galilean interactions allow for a successful 
screening of the force mediated by the Galileon on solar system scales, via the so-called Vainshtein mech- 
anism |8l], implies that these theories can in principle modify gravity in a way that is consistent with the 
most stringent experimental tests of GR. Nevertheless, the study of spherically symmetric backgrounds in 
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Galileon theories has revealed a number of undesirable features related to low strong coupling scales and 
superluminal propagation of fluctuations in the Galileon field. This has given rise to a wealth of studies 
that attempt to generalize the original Galileon model in a way that avoids these unpleasant issues while 

|<3l, DBI Galileon theories S, and multifield 



keeping its virtues. These include covariant completions 



generalizations 



HQ (see ta, 



1411 for reviews). 



lol . This model is the 



One particular model in this class is the bi-Galileon theory proposed in 
most general extension of the original Galileon theory to two scalar fields, in that it assumes only Galilean 
invariance (with second-order equations of motion) and a linear coupling between the Galileons and matter. 



16H analyzed the behavior of pertur- 



After developing some elegant algebraic methods, the authors of II15L 
bations on a spherically symmetric background with a massive point source, concluding that in this setup 
it is possible to avoid the issues of strong coupling and superluminality that affected the single-Galileon 
model. The goal of the present paper is to show that these conclusions need to be revised. In particular, we 
will show that superluminality cannot be avoided in bi-Galileon theory. We will also show that, even in the 
extension of this model to an arbitrary number of fields (what we will refer to as multi-Galileon theory), 
it is not possible to find a set of parameters that makes the theory free of superluminal propagation while 
avoiding the presence of extremely subluminal fluctuations and the related issues of strong coupling. 



When this work was in its final stages Ref. IU7I1 appeared. This paper has essentially the same scope as 
ours, and contains a proof of the presence of superluminal perturbations in the cubic multi-Galileon theory, 
as well as in the quartic bi-Galileon theory, in full agreement with our results. Our work can therefore be 
seen as a complement to that reference, since it studies the general multi-Galileon theory including both 



f. iHl 



cubic and quartic interactions. In addition, Ref. 01711 also studies the behavior of perturbations around a 
gas of particles with a given radial profile, finding that superluminal fluctuations cannot be avoided. This 
neatly illustrates that the phenomenon of superluminality is not a special feature of the configuration with a 
massive point source, but a very generic property of Galileon theories. 1 



II. MULTI-GALILEON THEORY 



The multi-Galileon theory considered here is perhaps the simplest generalization of the original Galileon 
model [1] to N fields. The theory is invariant under independent Galilean transformations in the fields, 
ir A — y 7Ta ~\~ b^Xf! + C A, with A = 1, . . . , N, and assumes a universal linear coupling between the Galileons 
and the trace of the stress-energy tensor. Then, by means of some trivial field redefinitions, we can always 
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condider a single Galileon, say 7ri, to be directly coupled to matter, so that the mixing between the Galileons 
and matter is encoded in the Lagrangian 



-7r, matter 



TTlT. (1) 



As it was done in 



... 



150, we make the working assumption that the contribution of the Galileons to the 
stress-energy tensor is negligible compared to the effects coming from /^matter- In other words, we neglect 
the gravitational backreaction of the Galileon fields. This allows for a consistent local analysis in which all 
fields can be taken to propagate in flat spacetime. For this theory to be an interesting and viable modification 
to GR, it is necessary to take into account the nonlinear interactions of the Galileons with themselves, 
allowing in principle for a successful Vainshtein screening on solar system scales. These interactions are 
encoded in the Galileon Lagrangian, given in 4 dimensions by 

5 

C"k — ^ ^ An (2) 

OQQ, 



where 



18] 



n=l 



£"n — ^ ] { a mi,...,m N 7T l + " " " + a mi,...,m N 7rN ) ^rni,...,m N j (3) 

m\-\ \-mN=n— 1 



and 



£ mi ,...,m N = (mi + • • • + m N )\S^ ■ ■ ■ daZ\ " " " 6 h " " " ^ 



(4) 



X [(a^TTi) • • • (a^^TTl)] • • • [{d Ul d^ N ) ■ ■ ■ (d„ mN d^NTT N )}. 

Explicitly we have 
£o,...,o = 1) 

£l,0,...,0 = Pi], 

«fi,i,o,..,o = Pi]p 2 ]-pin 2 ], 

fi,i,i,o,..,o = pi]p 2 ]p 3 ] - pi]p 2 n 3 ] - p 2 ]pin 3 ] - p 3 ]pin 2 ] + pin 2 n 3 ] + pin 3 n 2 ], 
i.:.„ „ = pi]p 2 ]p 3 ]p 4 ] - pi]p 2 ]p 3 n 4 ] - pi]p 3 ]p 2 n 4 ] - pi]p 4 ]p 2 n 3 ] 

- p 2 ]p 3 ]pin 4 ] - p 2 ]p 4 ]pin 3 ] - p 3 ]p 4 ]pin 2 ] + pi]p 2 n 3 n 4 ] 

+ pi]p 2 n 4 n 3 ] + p 2 ]pin 3 n 4 ] + p 2 ]pin 4 n 3 ] + p 3 ]pin 2 n 4 ] + p 3 ]pin 4 n 2 ] 
+ p 4 ]pin 2 n 3 ] + p 4 ]pin 3 n 2 ] + pin 2 ]p 3 n 4 ] + p 1 n 3 ]p 2 n 4 ] + pxn 4 ]p 2 n 3 ] 

- pin 2 n 3 n 4 ] - pin 2 n 4 n 3 ] - pin 3 n 2 n 4 ] - pin 3 n 4 n 2 ] - pin 4 n 2 n 3 ] 

- pin 4 n 3 n 2 ]. 

(5) 
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All the other Galilean invariants are obtained by exchanging or identifying different fields. We are using the 
notation (TIaYv = d^dv^A, and [II] denotes the trace of the matrix II. 

It is actually more convenient not to use the Lagrangian coefficients mN directly, but to use instead 

the coefficients appearing in the equations of motion as the parametrization of the theory. 2 These are given 
by 3 

d 



5ir A 



d xC n — (m-A + 1) I a mi ,...,m N + y ] a m 1 ,...,m B -l,---,mA+l,---,rn N 1 ^rni,...,rn N 

0<m.i+---+m JV <4 y B^A J 

V a A £ 

/ j u "mi,...,m N '-'rni,...,m N j 

0<miH \-m N <4 

(6) 

and we see that the Galileon coefficients mjv are related to the original Lagrangian coefficients by 



{m A + 1) [ adj,...,™^ + Yl a mi,-,m B -l,...,m A +l,..,m w ] ■ ( 7 ) 

B+A 



These coefficients are not all independent, however, since they satisfy some integrability conditions because 
of the fact that they are all derived from the same Lagrangian I15L The reader can check that, from Eq. ©, 
one obtains 

^fi"mi,...,mA- l,...,m.B,...,mjv ^Vi,...^,...^-!,...^' (^) 
The equations of motion for the Galileon fields can then be written as 

a m 1 ,...,m N ^m 1 ,...,m N = ~T, 

0<miH hmjv<4 

^ ] a mi,...,m N £mi,...,m N = 0, (A = 2, . . . , N). 

0<miH hmjv<4 



(9) 



In the following we will set the tadpole coefficients to zero. The reason for this is that we are 
implicitly thinking in field configurations on a self-accelerating de Sitter background; this was the initial 
motivation of the original Galileon model. The equations of motion for these configurations have precisely 
the same form as in ©, except that the constant terms do not appear because of the background equations 

110. 



2 Here and in the following we will use essentially the same notation as the one introduced in IH5LIK 

3 In the following expressions, ^ m is defined as zero whenever one of the indices equals — 1. 
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III. SPHERICALLY SYMMETRIC SOLUTIONS AND PERTURBATIONS 
A. Static spherically symmetric solutions 

For static, spherically symmetric field configurations we assume tta = t^a{t)- The Galilean invariants 



£ mi ,...,m N then reduce to 



£l,0,-,0 = -2-r( r TTl), 

2d. . ,. 
£i,i,o,...,o = -2-r{r'K l -K 2 ), 

r dr (10) 
2 d 



£l,l,l,0,-,0 = jT^^s). 



^1,1,1,1,0,. ..,o — 0, 

and all the other terms can be obtained by exchanging or identifying different fields. Note that the prime 
denotes differentiation with respect to r. We will focus on the case of a massive point source, for which 
T = — Mi5 3 (r), with M the mass of the source. The equations of motion © then become 



~(r 3 F\ yi ,...,y N ))=MSHr), 



r 2 dr 

~(r*F A (y u ...,y N ))=0, (A = 2,...,N) 



(ID 



where 

F A ( yi , ...,y N ) = f A + 2f A + 2f A , (12) 

fn = ^2 a m 1 ,...,m N y\ 1 ' ' ' Vn N ^ (13) 

miH \-m^f=n 

and we defined yA = K A /r. Integrating Eq. (fTTT t gives the following N algebraic equation for y±, . . . , y^: 

M 



F 1 (yi, ...,y N )- 
F A ( yi ,...,y N ) = 0, (A = 2,...,N). 



(14) 



Notice that both the coefficients mN and the variables yA are dimensionful. The scales of the Galileon 

coefficients are chosen to be such that, on the de Sitter background, the nonlinearities in the fields become 
important on Hubble scales and, in addition, that it is consistent to neglect the gravitational backreaction of 
the Galileons. Although this is not relevant for our present puiposes, let us simply mention that this choice 



of scales results in a Vainshtein radius of ||1 , 



16] 



M 



M£,H 



1/3 



rv ~ [ -tttifi ) • (15) 
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With this in mind, we can now perform some trivial rescalings and make the Galileon coefficients m 
dimensionless and (genetically) of 0(1), and to make the variables dimensionless as well. With a slight 
abuse of notation we will keep using the same symbols though. The equations of motion then become 

3 

(16) 



F A ( yi ,...,y N ) = 0, (A = 2,...,N). 

We see that the size of the variables yA is directly tied to the ratio ry jr in this notation. Next we define the 
matrices S n with entries 

(S n ) AB = J-f*. (17) 
OVA 

With the help of Eq. ([8]) it is easy to show that the matrices S n are symmetric. Notice that £2, S3, and S4 
depend on the yA, and therefore on r, but Si is a constant. In fact, the matrix Si is the same matrix that 
appears in the Galileon kinetic Lagrangian: 

£2 = ~ X^Ei WVa^ttb. (18) 

A,B 

We will therefore require that Si be strictly positive definite in order to avoid ghost instabilities. Using 
these definitions, the Jacobian matrix of the functions F A (yi, ... ,un) can be written as 

U = Si + 2S 2 + 2S 3 . (19) 

A continuous solution for the yA exists provided detU 7^ for all r > 0. At large distances from the source, 
r S> ry, the yA are small and the equations of motion are dominated by the linear functions ff. Morover, 
in this regime we have detU ~ detSi > 0, from where we see that detU > for all r > since detU 
cannot change sign. 

At large distances, Eq. ( fTBl implies that yA ~ (ry/r) 3 , from where it follows that tta ~ ry jr. This is 
of course the expected behavior for the regime where nonlinearities are not important. At short distances 



r 



<C ry, on the other hand, the equations of motion are dominated by the cubic functions f A , and we 



genetically expect to have yA ~ ry jr. It follows that tta ~ r/ry in this regime, which shows that the 
Galileon fields are indeed screened at distances smaller than the Vainshtein radius. 



B. Perturbations of the spherically symmetric background 

The next step in our analysis is to study the behavior of perturbations on the spherically symmetric 
background discussed above. We let tta — >■ ^A + 4>a, where 4>A{t } v) is a small fluctuation. To quadratic 
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order the Lagrangian for the fluctuations can be written as 

£ = ^d t $ • Kdt* - \d r <S> ■ Ud r <P - ^d n $ • Vdn$, (20) 

where <J? = (<p\, . . . , ^jv)> an d K, U, and V are N x N matrices. The matrix U is the same matrix that was 
denned in Eq. (fT9l ). In the above equation 8q denotes the angular part of the gradient operator in spherical 
coordinates. The matrices K and V can be most easily computed from the equations of motion for the 
fluctuations, with the result 

K = ( 1 + l r -f) ( S i + 3S 2 + 6S 3 + 6S 4 ) , 
V 2 dr) 

Notice that, to avoid ghost instabilities in the fluctuations, we must require that the matrix K be positive 



definite for all r > 0. It was shown in Ref . 1 3] tl 
strong coupling scale on the spherically symmetric background, precluding a consistent analysis in terms 
of the Galileon fields at short distances from the source. Explicitly, assuming a Vainshtein radius given by 
Eq. (fT5l) . one can show that the strong coupling scale in the presence of quintic Galileon interactions is 

A M~ rl/6 ^&i^ < 22 > 

Galileon excitations become strongly coupled when rA(r) ~ 1, which for the solar system corresponds 
to a critical distance r c ~ 10 5 km. Even though this distance is much smaller than the Vainshtein radius 
(ry ~ 10 16 km for the Sun), it is clear that the domain of applicability of the theory is quite limited, in the 
sense that the perturbative description in terms of the fields tta can no longer be trusted below the distance 
r c . We will therefore adopt the assumption that these quintic interactions are absent, so that £5 = 0. This 
also implies that £4 = in Eq. (|2TI ). 

The linearized equation of motion for the perturbations reads 

- Kdj§ + -^d r (r 2 Ud r $) + V8&$ = 0, (23) 

where 8q denotes de angular part of the Laplacian operator. For perturbations of sufficiently small scales, 
we can approximate these equations as 

- Kd?<S> + U8%$ + Vdl§ = 0. (24) 

In Fourier space we have 

[Ku 2 - Up 2 - Vpl] *(>,Pr,Pn) = 0, (25) 



8 



where p r and pq are, respectively, the momenta along the radial and orthoradial directions. Parametrizing 
the momenta as p r = pcosq, po, = psinq, we find that the squared sound speeds (? A {q) {A = 1, . . . , N) 
are given by the eigenvalues of the matrix 

M(q) =K~ 1 U cos 2 q + K^V sin 2 q. (26) 

Absence of gradient instabilities requires that c\(q) > for all r > and for all q. This means that the 
matrix M(q) must have nonnegative eigenvalues for all r and q. Absence of superluminal perturbations 
requires that c 2 A {q) < 1 for all r and q. This means that the matrix M(q) — I must have nonpositive 
eigenvalues for all r and q. 

One further requirement we would like to impose on the theory is that extremely subluminal modes 
with c 2 A <C 1 be absent. The reason for this is that, as shown in Ref. [1], this phenomenon of extreme 
subluminality is directly associated with a lack of damping in the fluctuations which results in a low strong 
coupling scale. This is undesirable at the theoretical level, since it significantly limits the domain of validity 
of the perturbative analysis; outside this domain the Galileon fields ir A cannot be reliably identified as the 
correct low energy degrees of freedom of the theoiy. Furthermore, there is also the phenomenological issue 
of the Cerenkov radiation associated with, for instance, the motion of the Earth in the field of the Sun. This 
has the implication that the Galileon field around the Earth cannot be consistently computed in the static 
approximation due to the retardation effects related to the slow propagation of the Galileon fluctuations. 

In the following subsections we will study the behavior of perturbations in the large- and short-distance 
limits. We will show that the above requirements cannot all be satisfied simultaneously. This is the main 
result of this paper. 



C. Behavior of perturbations at large distances 

We begin by studying the behavior of fluctuations at large distances from the source, distances much 
larger than the Vainshtein radius ry. This is the regime where the variables y A are small, and therefore 
the equations of motion are dominated by the linear terms f^. Following Jl6 1. we perform an asymptotic 
expansion in decreasing powers of r of the form y A = yy + y^ A + . . .. Let us assume first that neither the 
linear terms fi{y A ) vanish (this is equivalent to assuming that the kinetic terms for the fluctuations do not 
vanish), nor the quadratic terms f^(y A ) vanish. From the equations of motion we then find 

(0) 1 (1) 1 

Expanding the matrices S n perturbatively in a similar fashion, S n = £^ + + . . ., we find 

4 0) ocl, ... (28) 
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The matrices K, U, and V are approximately given by 

K ~ Si -3eJ, 1} -6S^ 0) , 



U ~ Si + 2S[, 0) + 2 (s^ + S^ 0) ) , (29) 



F-Sx-Sf -4(S«+S 3 



L > + E<^ 

where terms of order (ry/r) 9 were neglected. The matrix M(q) is given by 

M(q) ~ /+ (3cos 2 g- ljE^E^ + (6 cos 2 g - ljS^E^ + (6 cos 2 9 + 2)E^ 1 eJ )) . (30) 

We see that the matrix M(q) — I ~ (3 cos 2 g — l^E^Eg changes sign at cos 2 q = 1/3, implying that 
perturbations along some directions will be superluminal. We can avoid this by choosing coefficients such 
that E^ = 0. This in turn implies that the quadratic terms /^(y^) of the equations of motion vanish (see 
Appendix |A]), and so one must repeat the analysis starting with this assumption. 

Assume that the linear terms fi(y^ ) and the cubic terms /^(y? ) do not vanish, but that the quadratic 
terms f^(y^) as well as the matrix E^ do vanish. Then we find 

(°) 1 (!) 1 

»V K ^ va x ^ ••• ( 31 ) 

S^ocl E«ocl ... (32) 

The matrices K, U, and V are approximately given by 

K~ Ei - 6E^ 0) - 6E^, 

C/~E 1 +2Ef +2Ef } , (33) 
V r ^Ei-4E5° ) -7E^ 1) , 

where terms of order (ry/r) 12 were neglected. The matrix M(q) is given by 

M(g) ~ / + (6 cos 2 q + 2)E^ 1 T, { ° ) + (9 cos 2 q - ljE^E^. (34) 

It is clear that, for the matrix M(q) — I to have nonpositive eigenvalues, we need that the matrix Eg ^ be 
negative semidefinite. This requirement would make perturbations (slightly) subluminal at large distances, 
and in addition to the condition that Ei be positive definite, the theory would be free of instabilities in this 
regime. 
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D. Behavior of perturbations at short distances 

We repeat the analysis of the previous subsection, this time in the region of short distances from the 
source, r<rj/. Again we perform an asymptotic expansion, y& = + Va + • • •> tn i s ^ im& m increasing 
powers of r. We can then solve the equations of motion order by order, finding that 

y { A a v'a k1 i • • • ( 35 ) 

Note that this assumes that the cubic terms f^{y^ ), which dominate the equations of motion at short 
distances, as well as the quadratic terms f^(y^ ), do not vanish. Expanding the matrices E n perturbatively, 
i.e. E n = Si 0) + E^ + . . ., we find 

Sfocl, E^Efoci, ... (36) 
The matrices K, U, and V are approximately given by 

K ~ 2S( 0) + 2 (2E« + Ef ) , 

V~2Yf +2(Ef ) + S< 0) ), (37) 
F.EW+Ef, 

where terms of 0(1) were neglected. It is clear that K~ 1 V = 0{r /ry) at small distances, implying that 
there exist fluctuations in the orthoradial direction that are extremely subluminal when r <C ry (this is 
assuming the matrix K _1 V has only nonnegative eigenvalues; otherwise there is an instability). The way to 
avoid this would be to choose the parameters of the theory so that Eg ^ = 0. From the results of Appendix 
lAl this is not consistent with the assumption f% (y^ ) ^ 0, so we must repeat the analysis stalling from the 
assumption that both Eg ^ = and f^iy^) = at r < ry. 4 The equations of motion then imply that 
the da do not go as 1/r at short distances. Instead, assuming that the quadratic terms f% do not vanish 5 
(otherwise we would be left with the linear terms only, and so the Vainshtein mechanism would not work at 
all), we find 6 

yf^^ Va^^ ■■ (38) 

4 In principle it could be that f^{y2 ) — but £3 7^ 0. It is easy to see, however, that in that case the matrix K would go as 
l/r 3//2 at short distances, whereas the matrices U and V would go as 1/r 3 , implying the existence of squared sound speeds that 
are very large and positive or very large and negative, neither of which is desirable. 

5 In Appendix [D] we show that, for the case of two Galileons, the condition that £3°' — at large distances implies that the 
functions f A must vanish in the short-distance limit, implying that there can be no Vainshtein screening of the Galileons for 
r C ry- The following special case therefore applies to multi-Galileon theory with N > 3. 

6 Notice that to derive that oc 1 we used the fact that (/3*)' 1 ' = 0, as follows from the results of Appendix lAl 
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This gives 



1 



and the matrices K, U, and V are approximately given by 



[/^2fS^ 1) + 4 0) ) , (40) 



1 



J 3 ~ ^2 
^(1) , v (°) 



V~-[ E^ 

where terms of 0(1) were neglected. The matrix M(q) is then given by 

M a | cos* , (e« + Ie<»») (e« + Ef ) + i sin* , (E« + iEf (E« + E<°>) . (41) 
Consider the matrix 

A S (E« + isf)" 1 (E« + Ef) = 2/ - (E« + ^Ef)" 1 E« (42) 

From the results of Appendix |B] we have that detEg 1 = 0, from where it follows that the matrix A has 
at least one eigenvalue equal to 2. This implies that the matrix M has an eigenvalue equal to 4/3 when 
sin q = 0, corresponding to a superluminal mode along the radial direction. 

We conclude that by tuning the Galileon coefficients in such a way as to avoid the extremely subluminal 
perturbations that generically appear near the source, we end up in turn with a superluminal mode. It seems 
preferable, then, to go back to the assumption that f^(y^ ) ^ 0, which implies that the matrix E 3 °^ does 
not vanish identically, and simply accept the presence of an extremely subluminal mode (forgetting about 
the undesirable phenomenological and theoretical consequences of this). We still need to make sure, of 
course, that there are neither unstable nor superluminal modes. It is clear, from Eq. (l37T i. that we can 
avoid a gradient instablity and guarantee the existence of continuous solutions by choosing parameters such 
that Eg ^ is positive semidefmite. Notice that, by the results of Appendix we cannot have E 3 °^ strictly 
positive definite for r ry, since the absence of superluminality at large distances requires E 3 °^ to be 
negative semidefmite for r 3> ry. The only possible loophole would be that, by a very special choice of 
the Galileon parameters, Eg ^ results to be singular positive semidefmite (but nonvanishing) for r <C ry, 
and that it results to be singular negative semidefmite for r S> ry. We have no proof that, for an arbitrary 
number of Galileon fields, this special case also suffers from instablities or superluminality (in Appendix |E] 
it is shown that this loophole leads to a contradiction in the case of two Galileons). The main conclusion 
remains however: extremely subluminal perturbations will be present at distances close to the source. 
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IV. FINAL REMARKS 

Let us briefly summarize our results. We have studied the behavior of small fluctuations on static, 
spherically symmetric backgrounds in multi-Galileon theory. By imposing the requirements of stability, 
successful Vainshtein screening, absence of superluminality, and absence of extremely subluminal fluctu- 
ations, we have derived constraints in the large- and short-distance regimes. We have shown that these 
constraints cannot all be satisfied simultaneously in multi-Galileon theory and, in particular, that superlu- 
minal perturbations will always be present in bi-Galileon theory. It is an interesting open problem whether 
superluminality can in principle be avoided in multi-Galileon theory with N > 3, thanks to the loophole 
we identified at the end of Sec. Hill It seems that neither our results of Appendix [E] nor the proof given in 
Inl can be easily generalized to more than two Galileons. Perhaps the geometrical picture we developed 
in Appendix ICl will prove useful to solve this question. 

One may argue that the issues related to the retardation effects due to the extremely subluminal fluctua- 
tions are not serious, since after all the Galileon fields are screened at short distances and the corresponding 
Cerenkov radiation may not present a problem from the experimental point of view lltll . Of course, there is 
still the question of whether such a model can be consistently applied in the static approximation to the solar 
system, but this is not directly related to the consistency of the model at the theoretical level. The problem, 
however, is that the slow propagation of these fluctuations implies a lack of enhancement of the kinetic 
Lagrangian relative to the interactions, with the consequence of a very low strong-interaction scale for the 



fluctuations 



Perhaps more serious is the presence of superluminal perturbations, which were shown to 
be unavoidable in the bi-Galileon theoiy, and that will generically be present also in multi-Galileon theory. 
On the one hand, superluminality is a sign that the theory cannot be regarded as the effective low-energy 



description of a microscopic Lorentz-invariant theory 111 911 . However, this by itself does not mean that the 
theory is inconsistent, of course, unless one shows that closed timelike curves exist in the regime where 
the effective field theory is applicable. Interestingly, it has been conjectured that this does not happen 
in Galileon theory J2O], although it seems that further investigation is needed before we can reach a full 
understanding of the physical origin of superluminality in this class of theories. 
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Appendix A 

In this section we show that the condition £^ = implies that f A (y^) = for all ^4 = 1, . . . , N. 
From the definition 

fn = a m 1 ,...,m N y\ 1 ' ' ' Vn N ^ (Al) 



mi H \-m.N=n 



we have 



Therefore 



VB^-fn = £ ^Bai^^yT ■ ■ ■ VT- (A2) 

miH hmjv=n 



iV ^ 

^ZvB-g—fn = £ a m 1 ,...,m N yi ll ---yiV JV ( m l + --- + m Ar) 

_B=1 ^ B miH hm JV =n 

A mi . . . „, m N (A3) 



EA mi m 

a rni,...,m N Dl '"Vn 



n 

miH hmjv=w 



It follows that the condition 



(A4) 



(0) 



implies that f A {y^) = for all A = 1, . . . , N. Conversely, the assumption that f A (y^) f° r some A 
implies that the matrix £„ cannot vanish identically. 

Incidentally, notice that the condition £^ = also implies that the next-to-leading functions {f A )^ 
also vanish. Indeed, 



N 



N 

Y j y ( B ) ^ i n ) )BA = Q. (A5) 



y=y(0) B=l 



Appendix B 

In this section we prove that det = if S^ 0) = 0. Recall that we are working with an asymptotic 
expansion (either at large or short distances from the source) of the form i/a = + y^P + . . .. We can 
assume that all the y^ are proportional to each other, y^ = aAZ^°\ and likewise for the yj^ , that is 
yy = /3az^. The cannot be all equal to zero, and in particular we can set one of them equal to 1. The 
same holds for the Pa- Recall the definition 

f) 

/■^ \ -fA \ ■> A mi mn — 1 mjv / Dn 

(£„) Bj4 = — /„ = 2^ m sa mu ..., mjv yi • • -y B B ■■■y N - (Bl) 

miH hmjv=ra 
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We are assuming that E„ vanishes, and therefore 



(4 0) 



mi H \-rriN=n 



771 /v n 

■a^ = 0. 



The matrix £^ is given by 

= E 

miH hmjv=^ 



'mi,...,mjv 



(1) (0)' 



mi— 1 



(0)ms-l (0) mjv 
' "b " ' i/jV i" " 



. (1) (0)mi (0) 

+ m N yyy{' ■■■y B ' 
= (z®) n - 2 zW £ 



,v,m»-l (0)mjv — 1 



'mi,...,mjv 



miH hmjv=n 

+ (mg - IJPBai ■ ■ ■ a s " ' a TV ^ 1" m NPN<y-i 

Consider the linear combination 

TV TV 

E^(4 1) w = (. (0) r-v i )^ e 

B=l \-m N =n 

, n Wl 771m— 1 /o mi 

+ m A r/37 V a 1 • • • - Pbc4 • • • 

miH |-mjv=ra 



/O mi — 1 rriR — 1 m/vi 

mipic^ • • • a B H • • • a^y™ + 



mi mn— 1 mw — 1 

■ ■ ■ a- • • • ™.. iV 



B 



a 



TV 



m B a mi ,...,m N 



a mi— 1 mw i 



mj-l mjv 
a B ' ' ' 01 TV 



.,mjv 



a mi— 1 mjv , 

m\P\a x ■■■a N " + ■ ■ 



i o mi m/v- 

+ m A rp/va 1 1 • • • a/ 

TV 

-(z^r-w^ft* E 

B=l miH hmjv = 



m BV,,miv l Q 



1 ' B ' ^ 



TV 



(B2) 



(B3) 



(B4) 



= 0, 

from Eq. (IB2I) . This shows that the rows of the matrix are not linearly independent, and therefore 



det = 0. 



Appendix C 



In this section we prove that the matrix S3 (thought as a function of the N variables yX) cannot be 
positive definite at one point, and be negative semidefinite at some other point. The key observation is that 

fi 



the can all be derived from a single function L3. To see this, start from 



dL 3 
dyi 



fl= E 



mi „, m N 



l mi,...,m N yi ' ' ' Vn ' 



(CI) 



miH hmjv=3 
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and integrate to find 

a 1 

j \ ^ u 'mi,...,mjv mi+1 m N . r , \ 

Ls= E mi + 1 Vi ~'Vn +9i{vA^yx) 

miH hmjv=3 

a 1 , ^ C2 ) 

E^mi — l,...,mjv mi m N , / / \ 
2/1 ■■■y N +9i{va + yi)- 
mi 

miH hmjv=4 

Here 51 is some function that depends on all the ua except y\. Differentiating with respect to y 2 we obtain 

9L 3 \ - a ln 1 -l,...,m N _ mi m 2 -l „mjv , ^1 



dy 2 A m i ^2/2 

miH hmjv=4 

mi^O 



mi H h"ijv=3 

mi^O 



E2 mi mjv i 

a mi,...,m N yi Vn + Qy 2 ' 



(C3) 



where we used Eq. dU). Equating this with /f , we get 



5y 9 — X/ ^^....mjvyr 1 " ■ Vn n ■> 



(C5) 



miH hm.jv=3 

mi=0 

which is independent of y\. Integrating again gives 

a 2 

E^mi,...,m N mi m 2 +l mjv i / _L \ 
rrVi % ■■■fiv + 92{yA^ yi,y2) 
m2 + 1 

miH hmjv=3 

mi=0 

a 2 

\ r m\,m2 — l...,mjv mi mw . / / \ 

miH hmAr=4 

mi=0,m2^0 

where <?2 is a function that depends on all the j/a except j/i and 1/2- Repeating this process N times finally 
gives the function L 3 (up to an irrelevant integration constant): 

n 1 n 2 
T _ \ ~* u m 1 -l,...,m N mi m N , \ ~» u, m 1 ,m 2 -l...,m N mi „ mjy 

-^3 — 7> y i '"Vn + 7> y i " " " 2/jv 

^-^ mi ' mo 

miH hmjv=4 miH \-mtq=A 



mi^O mi=0,m 2 ^0 

AT (Co) 

+ ...+ £ vr • • • //v x ■ 

m jv 

miH hmjv=4 

mi=0,...,mjv-i=0,mjv7^0 

From Eq. © we see that the coefficients appearing in these sums are all independent, and so the complete 
sum can be written as 



L, 



E 

miH hm JV = 



.4 



mi , 



.,m N y± 



UN ' 



(C7) 
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where the A mi mN are all independent parameters. Thus, the function L3 is a general iV-ary quartic form 
(a homogeneous polynomial of degree 4 in N variables). The entries of the matrix S3 can then be written 

as 

(S 3 )ab = ^-fi = t^~ L ^ ( C§ ) 
ay a oyAOyB 

that is, S3 is given by the Hessian matrix of L3. Thus, our task is to prove that the Hessian matrix of a 
general quartic form cannot be positive definite at some point and negative semidefmite at some other point. 

We first show that this is true for N = 2 (the case N = 1 is trivially true). We have to show that the 
Hessian matrix of a general binary quartic form, 

q(x, y) = ax A + bx 3 y + cx 2 y 2 + dxy 3 + ey A , (C9) 

cannot be positive definite at one point, say (xi,yi), and negative semidefmite at some other point, say 
(^2,2/2)- The proof of this is greatly simplified by the fact that any quartic form can be reduced to its 



canonical form (see for example Il2lll ). 

Q(X,Y) = rX A + §mX 2 Y 2 + sY A , (CIO) 

by means of a nonsingular linear transformation (x, y) 1— s> (X, Y). The Hessian matrix of Q(X, Y) is given 
by 

/ rX 2 +mY 2 2mXY \ 
H Q (X,Y) = 12 \ (Cll) 

V 2mXY sY 2 + mX 2 J 

Assume that there is a point (X\ , Y\) ^ (0, 0) where Hq is positive definite. 7 Then its eigenvalues must be 
strictly positive, or equivalently, its trace and determinant must be positive: 

(rX 2 + mY 2 ) + (sY 2 + mX 2 ) > 0, (C12) 

(rX\ + mYl){sYl + mXf) - 4m 2 X 2 Y? > 0. (C13) 
In particular, these imply that 

(rX\ + mY?) > 0, (sY? + mXf) > 0. (C14) 

Assume also that there is another point (^2,12) 7^ (0,0) where Hq is negative semidefmite. Then its 
eigenvalues must be nonpositive, or equivalently, its trace must be nonpositive and its determinant must be 
nonnegative: 

(rX 2 + mY?) + ( S y 2 2 + mXl ) < 0, (C15) 



7 The condition (Xi, Yi) 7^ (0, 0) comes from the requirement that the background Galileon fields do not both vanish. 
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[rX\ + m,y 2 2 )(sY 2 2 + mJff ) - 4m 2 X 2 Y 2 2 > 0. (C16) 
In particular, these imply that 

(rX\ + mF 2 2 ) < 0, (sY 2 2 + mXf ) < 0. (C17) 

It suffices to consider the cases where (r = 1, s = 1), (r = — 1, s = 1), and (r = 0, s = 1) (we leave to 
the reader to check that in the cases where both r and s are zero, or when m = 0, one easily arrives at a 
contradiction). 

If r = 1, s = 1, then Eqs. (ICl~2l and (ICBT ) give 

(l + m)(X 2 + >f) >0, (l + m)(X 2 2 + y 2 2 ) <0, (C18) 

which is a contradiction. 

If r = -1, s = 1, then Eqs. (|CT4l > and dUTTT ) give 

X 2 

- X 2 + mY? >0 m>— i->0, (C19) 

l 

y 2 

y 2 2 + mXf <0 ^ m < ^^0, (C20) 

X 2 

and we arrive at a contradiction. 

If r = 0, s = 1, then Eqs. (IC141) and (ICTTT) give 

raFi 2 > 0, mlf < 0, (C21) 

which implies m > 0, y 2 = 0. But if that is the case then Eq. (IC 15b gives 

mXf < 0, (C22) 

and we arrive again at a contradiction, since X 2 and y 2 cannot be both zero. 

We conclude that the Hessian Hq cannot be positive definite at one point and negative semidefinite at 
some other point. Since positive (or negative) definiteness of the Hessian is a statement about the local con- 
vexity (or concavity) of a function, and since convexity is preserved by linear transformations, we conclude 
that the Hessian of the general quartic form q(x, y) also has this property. 

Finally we prove the general case of N variables using induction, by exploiting the relation between 
the positive definiteness of the Hessian and convexity. Assume that a general iV-ary quartic form cannot 
be concave at one point if it is (strictly) convex at some other point. But assume that there is an (N + 1)- 
ary quartic form Q that is (strictly) convex at a point x' = (x[, . . . , x' N+1 ), and concave at a point x" = 
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(x'{, . . . , Consider a hypeiplane containing the points x' , x" and the origin (such a hyperplane 

always exists for N > 2; this is why the case of one Galileon does not work as a base case for the inductive 
argument). Write the equation of this hyperplane by solving for one variable, say xa, in terms of the others 
(this equation is homogeneous, since the hypeiplane contains the origin), and consider the iV-ary quartic 
form Q' obtained by constraining Q to this hyperplane. By assumption Q' cannot be strictly convex at x' 
and concave at x" . But then Q cannot satisfy this property either. The reason for this is that a function / is 
convex (concave) in a region if and only if the function obtained from constraining / to any line contained 
in that region is convex (concave) also. In particular, for a function / of 3 or more variables, the function 
obtained from constraining / to any plane (contained in the region where / is convex) must be convex as 
well. We conclude that the (N + l)-ary quartic form Q cannot be strictly convex at x' and concave at x". 
This completes the inductive argument. 



Appendix D 

In this section we show that, in bi-Galileon theory, the special case in which the matrix £3 vanishes in 
the short-distance limit implies the absence of a successful Vainshtein screening of the Galileons. We will 
use the following notation for the Galileon coefficients: aj ni m2 = a mijm2 and a^ ll m2 = b mi:in2 . Notice 
that, with the exception of 601, 602 > and 603 > all the coefficients f> TOl , m2 can be expressed in terms of the 
coefficients a mi>m2 by using Eq. ©. 

Recall that we are working with asymptotic expansions at short distances, tja = + + ■ ■ ■> an d at 
large distances, yA = + V^A + • • •> where A = 1, 2 in bi-Galileon theory (here and in the following, the 
bar labels quantities evaluated in the large-distance limit, to avoid confusion with the short-distance limit). 
Recall also from Sec. [Ill] that we imposed the requirement that S^ ' = £2(2/1 ,2/2 ) = in order to avoid 
superluminal propagation at large distances. In this regime the equations of motion are dominated by the 
linear functions fi(y[ , y\ ), and the solutions are easily found to be 

d°) = ^ ( rV \\ y(0) = ( rV )\ (Dl) 

1 (aio&oi - all) ^ r ' 2 (010601 - «oi) ^ r ' 
The matrix i s given by 

■ (0) _ / 2a 2 o/3 + an aiifi + 2ao2 \ _( ) 

an/5 + 2a 02 2a 02 /3 + 26 02 

where we defined j3 = §1/^2 = ~ 601/001 (we will come back to the case where = later). The 
condition £3 = then implies that the cubic Galileon coefficients can be written in terms of 020 and j3 as 
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follows: 

a n = -2a 20 f3, a 02 = a 20 /3 2 , b 02 = -a 20 (3 3 . (D3) 

We are also assuming that, at large distances, the matrix £3°^ = S3 (yj°\ y 2 ) vanishes, as required to 
avoid extremely subluminal fluctuations. From the results of Appendix |A] we know that the cubic functions 
fzhn-ilh) a ^ so vanish, and so the equations of motion in the short-distance limit are dominated by the 
quadratic functions f^iVi ,y^)- Explicitly we have (see Eq. (fT6l) ) 

3 

J 



2f 2 \y?\y^) = 2 (a 20 « 2 + alia + a 02 ) yf 2 = (^Y 



2f 2 (y[ \yW) = 2 (°fc? + 2a 02 a + b 02 ) yf 2 = 0, 



(D4) 



and we defined a = yf ■* / yp (we assume for now that j4 0) / 0; we will come back to the case where 
y^ = later). However, if we use relations (|D3t we find that 

fl{yf\y < i ) ) = a w {a-pfy^\ 

fi(y? ) JV) = -a 20 H«-P) 2 yr- (D5) 
Since /3 7^ (from the requirement of having a positive definite kinetic Lagrangian for the perturbations), 
the second equation of motion in (ID4b implies that either a 2 o = (in which case all the cubic Galileon 



coefficients vanish) or a = (3. In either case we will have that /| (yf^, y^ ) = identically. This implies 
that the equations of motion will not be dominated by the quadratic functions f 2 {y^ , y 2 )> but instead they 
will be dominated by the linear functions ff(y^ , y^), meaning that there will be no Vainshtein screening 
of the Galileons for r < ry. 

The same conclusions arise in the cases where either y£ = or yij = 0. If yf = and y^ ^ then 
the short-distance equations of motion read 



(0)2 n 



(D6) 



which imply that an = 0. But then the condition E)^ = implies that all the cubic Galileon coefficients 
vanish and there can be no Vainshtein mechanism. If both y^ = and y%' = then this condition implies 
that a 2 o = and 002 = (with 602 unconstrained), and again there can be no Vainshtein mechanism. 
Finally, if y% 7^ and y^ = then the condition £3 = implies that 020 = 0, an = and 
a 02 = (while the short-distance equation of motion will require 602 = 0), and once again there will be no 
Vainshtein screening of the Galileons. 
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Appendix E 

In this section we show that, in bi-Galileon theory, the loophole mentioned at the end of Sec. ITTT1 leads 
to a contradiction after using the equations of motion. In the following we will use the same notation 
introduced in Appendix |D] At large distances the equations of motion are dominated by the linear functions 
fi~(Vx , y 2 )> an d the solutions are given by Eq. (ID 1 b . At short distances, on the other hand, the equations 
of motion are dominated by the cubic functions /^(y^ , y 2 °^)- Explicitly we have 



2/3^i°U 0) ) = (^)' 
2/3 2 ( y i°U 0) ) = 0. 



(El) 



The assumptions of the loophole case are that £g°^ = S3 (y ^ , y^ ) is singular positive semidefmite (but 
nonvanishing), and that £g°^ = £ 3 (y[°\y2°^) * s singular negative semidefmite. For the matrix £ 3 °^ to be 
singular we have two options; one is that the rows of this matrix are nonzero but proportional to each other, 
and the other is that one of the rows is zero. The first option does not work because of the equations of 
motion, for if the second row is A times the first row (with A 7^ 0), then 

p2^,(o) „(0)n _ 1 L(0), v (0K , fowo). 



hVUi 1V1 )- 3 Vi (^3 J2i+y 2 \ h 3 m 



A 
3 



V\ O3 hi + y 2 {^3 A2 



(E2) 



= A/ 3 i (yr,y 2 u 

where we have used the results of Appendix |A] We see that /f (yf^, y 2 °^) = (as implied by the second 
equation of motion) if and only if f^(y[ , y 2 ) = also, contradicting the first equation of motion. The 
only possibility is to have the second row of the matrix Eg ' 1 equal to zero (the first row cannot be zero, 
again by the equation of motion and the assumption that £3 is nonvanishing). With this choice the second 
equation of motion is satisfied as an identity. 

To analyze this case, it is convenient to work with the ratios a = yf^ /y^ and /3 = yf^ /y^ (this 
assumes 7^ and y 2 7^ 0; we will come back to this assumption later). The matrices and Sg ^ 
are then explicitly given by 

I 3a 30 a 2 + 2a 2 ia + ai 2 a 21 a 2 + 2a 12 a + 3a 3 1 

S3 = I I y 2 , (E3) 

a 2 \o? + 2ai 2 a + 3ao3 a\ 2 a 2 + Qa^a + 36q3 



- (0) _ / 3a 30 /3 2 + 2a 2 i/3 + a X2 a 21 (3 2 + 2a 12 f3 + 3a 3 \ _ (0 ) 2 

^3 — I I Vi ■ ( E4 ) 

a 2 i/r + 2a 12 (3 + 3a 03 a 12 f3 z + 6a 03 P + 3b 03 
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Sg ^ is singular positive semidefinite, with vanishing second row, if 



3a3oa + 2a 2 ia + a\ 2 > 0, 
a 2 ia + 2aua + 3ao3 = 0, 



(E5) 



a± 2 a + 6003a + 36| 



03 



0. 



(E6) 



£3°^ is singular negative semidefinite if 

(3a 30 /3 2 + 2a 2 i/3 + a 12 )(a 12 /3 2 + 6a 03 /3 + 36 03 ) - {a 2l p 2 + 2a 12 /3 + 3a 03 ) 2 = 0, 

(3a 30 /3 2 + 2a 21 /3 + a 12 ) + (ai 2 /3 2 + 6a 03 /3 + 36 03 ) < 0. 

Notice that we must obviously have p 7^ a to satisfy these conditions. Also, notice that /3 ^ 0; this follows 
from the explicit expressions we have for the large-distance and (see Eqs. (ID II) ). from where it 
follows that /3 = — boi/aoi, which is nonzero since 601 7^ from the condition that the kinetic Lagrangian 
be strictly positive definite. We will next show that conditions (IE5b and (IE6b are inconsistent. 

Suppose first that {a\ 2 p 2 + 6a 03 /3 + 3603) is strictly negative (Eq. (|E6b implies that it is nonpositive). 
Then we can solve for 3030 from the first equation in (IE6K finding that 



3a 30 = — 



P 2 



In addition, from Eqs. (IE5I ) we have 



(a 21 /3 2 + 2a 12 (3 + 3fl 03 ) 2 
(ai2/3 2 + Qa 03 (3 + 36 3) 



(2a 2 i/3 + ai 2 ) 



(E7) 



3ao3 = — «2i" — 2ai2a, 
3&03 = -ai2« 2 - 6a 3a, 

from where we find that 

(ai2^ 2 + 6a 03 /3 + 36 3) = (P - a)(a 12 (P + a) + 6003), 
(a 21 /3 2 + 2a 12 p + 3a 03 ) = ((3- a)(a 21 (/3 + a) + 2ai 2 ). 
Using Eqs. (lE7b and (IE9I ). and performing some straightforward manipulations, we find 



a 



(3a 30 a + 2a 21 a + a l2 ) = 



(P - a)(a 2 \{P + a) + 2a 



12 



(ai2(/3 + a) + 6003) 



(2a 2 i/3 + a 12 ) 



(E8) 



(E9) 



(2a 2 ia + CL12) 



(P-<x) 

P 2 
(P ~ «) 3 



a 2 {a 21 (P + a) + 2a 12 f 
{a 12 (P + a) + 6a 03 ) 
(a 21 a + a 12 ) 2 



+ (2a 21 aP + a 12 (P + a)) 



P 2 (a 12 (P + a) + 6a 03 ) 
(P - a) 4 (a 2 ia + ai 2 ) 2 
/3 2 (ai 2 /3 2 + 6a 03 /3 + 36 03 ) ' 



(E10) 
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where in the last line we used again Eq. (IE9I ). Since we have assumed that {a\ 2 f3 2 + 6003/? + 3&03 ) is strictly 
negative, we conclude that (3a3oa 2 + 2a2ia + a\ 2 ) must also be strictly negative, contradicting condition 
(ED. 

Let us next assume that (ai 2 p 2 + 6003/? + 3&03) = 0. Conditions (IE6I) then imply that also (a2i/3 2 + 
2ai2/3 + 3ao3) = 0, and that (3a^/3 2 + 2a 2 \f3 + 012) < 0. Taking these conditions together with (IE5I ). we 
see that a and f3 correspond to the two distinct roots of the polynomials 

Pi(x) = a 2 \x 2 + 2ai 2 x + 3a 3 = a 2 i(x — a)(x — 0), 

(Ell) 

P2(x) = ai 2 x 2 + 6a 3X + 36 3 = a\ 2 {x - a)(x - (3). 
(This assumes implicitly that a\ 2 7^ and a 2 \ 7^ 0, for otherwise we cannot have a ^ (3 and satisfy the 
above conditions.) Since the two roots of the polynomials are the same, we can find the following relations 
between their coefficients: 

a\ 2 = 3a 2 ia 3, 012003 = 021^03- (El 2) 

But if we now compute the discriminant of p\, we find that it is given by Ai = 4a 2 2 — 12a2iao3 = from 
the above relations. One can similarly show that A2, the discriminant of p 2 , also vanishes. This contradicts 
the fact that a ^ (3. 

Finally, let us drop the assumptions we made at the beginning about yy and y^ being nonzero. First, 
if y^ = 0, then the matrix E 3 °^ simplifies to 

S ( 0)= (0)2 (E13) 

\ a 2i ai2 J 

and we have to set 021 = = 012 and require 030 > to have E 3 ^ singular positive semidefinite and 
satisfy the equations of motion. The matrix E 3 °^ is then given by 

g (0) = / 3a 3 o y f )2 3a 032 /f 2 \ 

\ oa 3y 2 ba 03 y l y 2 + 6b 03 y 2 J 

For E 3 °^ to be negative semidefinite, we need Sa^y^ 2 < 0, which contradicts the condition 030 > unless 
y± = 0. But if this is the case then the equations of motion in the large-distance limit would imply that 



601 = 0; see Eq. (ID lb . This would contradict the requirement of having a strictly positive definite kinetic 
Lagrangian for the perturbations. 

Next, let us assume y 2 ^ but y 2 = 0. The matrix E 3 °^ is then given by 

S (0) = ( 3 ° 30 ° 21 ) y f )2 . (E15) 
a 2 i a 12 
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The condition that this matrix must be singular negative semidefinite implies 

3a3oai2 = O21, 

a 30 < 0, (E16) 
012 < 0. 

Since ^ 0) / 0, as we assumed above, the Galileon coefficients must satisfy conditions (1E5I) . If a\i < 0, we 
can solve for 030 from the first of Eqs. (IE16K finding that 

n 2 

2 21 2 

3a3oa + 2a 2 ia + ayi = a + 2a2i« + 012 

012 (E17) 

= (a 2 ia + ai2) 2 < 0, 

ai2 

contradicting conditions (IE5b - If ai 2 = 0, then from (IE16b . 021 = also. The first equation in (IE5b then 
reduces to 3a^a 2 > 0. This in turn implies that 030 > 0, contradicting (IE16I ). 
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